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We consider a nematic liquid crystal confined by two parallel flat substrates whose anchoring 
conditions vary periodically in one lateral direction. Within the Gaussian approximation, we study 
the effective forces between the patterned substrates induced by the thermal fluctuations of the 
nematic director. The shear force oscillates as function of the lateral shift between the patterns on 
the lower and the upper substrates. We compare the strength of this fluctuation-induced lateral 
force with the lateral van der Waals force arising from chemically structured adsorbed monolayers. 
The fluctuation-induced force in normal direction is either repulsive or attractive, depending on the 
model parameters. 

PACS numbers: 61.30.Dk, 61.30.Hn 
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I. INTRODUCTION 



Liquid crystals are characterized by large thermal fluc- 
tuations in their local orientational order arising from 
collective alignment of the long axis of their constituent 
molecules Due to such soft anisotropy, liquid crys- 
tals tend to respond easily to external forces. Confin- 
ing geometries such as thin films, on which most ap- 
plications of liquid crystals are based, change the fluc- 
tuation spectrum. This can cause not only structural 
changes 0, Q but also leads to fluctuation-induced ef- 
fective forces between the substrates Q, also known as 
thermodynamic Casimir effect. In correlated fluids such 
as liquid crystals this Casimir force exhibits a universal 
power-law decayas a function of the separation between 
the substrates jfj. However, this behavior is modified 
in the presence of other characteristic scales in the sys- 
tem 0,0, • In the case that the substrates are laterally 
modulated, discrete lateral modes of thermal fluctuations 
are also excited. Under such conditions, in addition to 
the forces acting perpendicularly to the substrates, effec- 
tive lateral forces arise 0, ^| with potentially interest- 
ing technological applications. We study the influence of 
anchoring conditions, which vary periodically in one lat- 
eral direction, on the fluctuations of a uniformly ordered 
nematic liquid crystal. Obviously, the periodicity £ of 
the substrate pattern gives rise to an oscillatory behav- 
ior for the lateral force as a function of the lateral shift 
5 between the substrates. For small inhomogeneities the 



'Present address: Institute for Studies in Theoretical Physics 
and Mathematics (IPM), School of Physics, PO Box 19395-5531, 
Tehran, Iran 



lateral force is proportional to sin(27r<5/£). (The analysis 
of nonperiodic patterns would provide an understanding 
of nematic phases exposed to chemically disordered sub- 
strates.) The present study actually extends our previ- 
ous work p"H where we considered the case in which only 
one of two confining substrates exhibits a chemical pat- 
tern so that there are no lateral forces. Here, in addition 
to the fluctuation-induced lateral forces, we calculate the 
lateral force between the patterned substrates across the 
vacuum, i.e., the background van der Waals force acting 
parallel to the substrates. This background force is gen- 
erated by the necessary chemical modulations providing 
the laterally varying anchoring strengths. 

In Sec.[H]our model and the theoretical formalism are 
specified. In Sec. IIII Al the fluctuation-induced lateral 
force is obtained. In Sec. IIII Bl we calculate the lateral 
van der Waals force between the patterned substrates. 
The results for the fluctuation-induced normal force are 
presented in Sec. IIVI and finally Sec. summarizes our 
results. 



II. SYSTEM AND FORMALISM 

We consider a nematic liquid crystal confined by two 
flat but chemically patterned substrates at a separation 
d. The patterns on the both substrates consist of the 
same periodic stripes of anchoring energies per area W a 
and Wb along the x direction but shifted relative to each 
other by the length 8 (see Fig. 1). The substrates are 
translationally invariant in the y direction. The stripes 
are considered to vary with respect to the strength of 
homeotropic anchoring so that the mean orientation of 
the director n is spatially homogeneous but the thermal 
fluctuations vary laterally giving rise to effective lateral 



2 




W a 



FIG. 1: The geometry of the nematic cell with patterned 
substrates. The patterns on both substrates are the same 
but shifted relative to each other. The patterns consist of 
periodic stripes of anchoring energies per area W a and Wb 
with the widths C, a and C,b, respectively. The wavelength of 
the periodicity is denoted as ( — £ a + £ f, and the lateral shift 
between the origins of the patterns on the top and the bottom 
substrate is denoted by S. Anchoring at both boundaries is 
homeotropic everywhere so that the thermal average of the 
director field no = z is spatially homogeneous. 

forces. Based on the bulk structural Frank free energy Q 
given by 



F h \n] = - d 



Ki{V ■ n) 2 + K 2 (n- V x n) 5 
+A' 3 (n x V x n) s 



(1) 



where V is the nematic volume, K\, K 2 , and K$ are 
the splay, the twist, and the bend elastic constants, re- 
spectively, the free energy of Gaussian fluctuations in the 
one-constant approximation reads 



W b [l — a(x)] at the lower substrate located at z = and 
W z=d (x) = W a b(x) + W b [l-b{x)} at the upper substrate 
located at z = d where 



a(x)= J2 e(z-fcc + £)e(fcc + 



and 



b(x) = a(x + S) 



(4) 



(5) 



describe the stripe modulations at the lower and the up- 
per substrates, respectively, <d(x) is the Heaviside step 
function, and £ is the periodicity. The stripes have the 
same width £ a = ( b = C/2 and are separated by sharp 
chemical steps. The functions a(x) and b(x) equal one 
at the regions characterized by W a and zero elsewhere at 
the lower and the upper substrates, respectively. Thus 
the surface free energy of the Gaussian fluctuations given 
by F s [ui,V2] = F a [v{[ + F s [v 2 ] reads 



Fr>] = \ Wa I d 2 x [i/(x, z = 0)] \(x) 



+W b / d 2 x [i/(x, z = 0)] [1 - a(x)} 
Js 



(6) 



at the lower substrate and 
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z—d r 



W a / d 2 x[v(*,Z = d)] 2 b(x) 



s 



+W b / d 2 x [i/(x, z = d)] [1 - b(x)\ 
Js 



(7) 



at the upper substrate. 

Minimization of total free energy F = -Ft, [^1,^2] + 
F s z=0 [^ 1 ,i/ 2 ] + F*~ d [v\, v^] leads to two boundary con- 
ditions: 

- Kd z y(x, z) + W a v{~x-, z) a(x) 

+ W b is(x, z) [1 - a(x)] = 0, z = 0, (8 a) 



Fb[vi,V2] 



t- 

^ 2 



d 3 a;[V^(x,z)] 2 , 



(2) 



where i = 1,2, is either of the two independent com- 
ponents of the fluctuating part 5v — n— no of the director 
n, K = K\ = K 2 = K3 is the effective elastic constant, 
and x = (x, y) are the lateral components of the Carte- 
sian coordinates r = (x, z). 

To describe the interaction of the liquid crystal and the 
substrates, we employ the Rapini-Papoular surface free 
energy given by 



(3) 



F a [n] = ~ f d 2 x^ 2=0 (x)(n-z) 2 - 
2 Js 

1 [ A 2 xW z = d {^){n-if =Fr° + Fr d 

2 Js 



where W is the anchoring energy per area S and z is the 
unit vector in z-direction. Here W z=0 (x) = W a a(x) + 



Kd z v(-K, z) + W z/(x, z) b(x) 

+ W b v{x, z) [1 - 6(x)] = 0, z = d, (8 b) 



where v is either v\ or 1^2 ■ 



A. Path integral technique 

The normalized [see, c.f., after Eq. JJJJl] partition func- 
tion Z of the fluctuating fields v^, i = 1,2, subject to 
the boundary conditions given by Eqs. (|8 aj l and <|8 bl) 
can be calculated within the path integral approach (see 
Ref. and references therein). The boundary condi- 
tions act as constraints which can be implemented by 
delta functions. They, in turn, can be written as inte- 
gral representations by introducing two auxiliary fields 
localized at z — and z — d, respectively [T2I fl3| . Af- 
ter performing the corresponding Gaussian integrals over 
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i>i, i = 1,2, the path integral reduces to a Gaussian 
functional integral over the auxiliary fields with a ma- 
trix kernel M, so that obtaining the result for Z reduces 



J 



to calculating (det M) -1 / 2 . For the geometry considered 
here, the matrix M is found to have the following matrix 
elements M a ^, a,(3 = l, 2: 



M n 



M- 
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(x,x') = { 
Ab(A{, — A a ) 



1 H a(x) 



-, Aft — A a , 
1 H r a(x ) 

Aa 



i(x) - a(a/)]a» - A 2 <9 2 }G(x - x', 



Mia (x,x') 



Ab — A 

A Q 



[a(x') + b{x)] - 2A b d z 



Ab(Ab — A a ) 



A Q 



M 21 (x,x') 



A? ' A -) 2 a(a;')%)+A6a 2 ,) G(x - x'.z - z') 



2 — d.z'— 



1 + 



Ab - A 



-[a(a;) + &(a;')]-2A 6 & 



A&(Ab — A Q ) 



Ab- A a \ 2 



(x,x') = { 



Aa A a 

a{x)b{x') + A 2 <9 2 )] G(x -xf,z- z') 

Ab — A a , 



a(x') + &(#)R 



a(s) + &(x')R 



d,z— 



A n 



-I Ab — A a , , 

1 + — b(x ) 



A b (A b K) [b(xf) b{x)]dz _ X 2 d 2\ G(x _ x ^ z _ 7j) 

Aa J 



(9) 



r 



tion length and G(r,r') 
the two-point correlation 
Vi in the bulk with its 



where A a (b) = K/W a g,) is the so-called extrapola- 

= k B T/{4nK\r - r'|) is 
function of the scalar field 
statistical weight given by 
j y d 3 a;j/ i (x, z)V 2 ^(x, z)} where the nor- 
malizing factor Zq is the bulk partition function and k B T 
is the thermal energy. 

In terms of the partition function Z, the free energy is 
given by F = -k B Th\Z = (fcsT/2) In det M, We note 
that normalizing Z by Zq amounts to subtracting the 
bulk free energy, so that the free energy F includes only 
the surface free energy and the finite-size contribution. 
The surface free energy depends neither on d nor on 8, 
so the fluctuation-induced force T = —dF reads 



T =- 



k B T 



Tr (M~ 1 dM), 



(10) 



where d is either d$ or dd corresponding to lateral or nor- 
mal displacements giving rise to lateral or normal forces, 
respectively. 



B. Periodic modulation 

The matrix kernel M is a functional of the pat- 
terning function a(x) on the substrates and therefore 
calculation of the inverse of M is nontrivial. How- 
ever, in systems with in-plane symmetries one may pro- 
ceed by a lateral Fourier transformation with respect 



to the lateral coordinates x = (x,y). In Ref. |14j . 
it is shown how the electrodynamic Casimir force can 
be calculated for a periodically modulated substrate 
(see also Refs. 0, In this reference, the lat- 

eral periodicity is used to transform the matrix M to 
a block-diagonal form in Fourier space (p, q) in which 
M(p,q) = / / d 2 x dVM(x,xV px e iqx '. Similarly, 
also here the matrix elements of the block AL 



with Mj t ki(jp y ,q y ) = 2ir5(py + q y ) B 
1, • • • , N = L/C,, are given by 



V L ■ 



Py)> 3 = 



k-l 



27TJ 



= N n 



-k-l 



2ttZ 

T 



(11) 



for I, k S Z, where L is the lateral extension of the sys- 
tem in the x direction and the N m are (2x2) matrices 
providing the following decomposition of the matrix M: 



M(p,q) = (2tt) 2 5{p y + q y ) ^ N m { PxiPy ) 

m— — oo 

r / 27rmA 
X0\p x + q x + — 1 . (12) 



However, it is interesting to note that the matrix 
M(x, x ) = M(x — x , a(x), a(x )) can also be repre- 
sented, using a more direct derivation than in Ref. [l4| . in 
a form in which M is diagonal fl6| . In view of the discrete 
lateral periodicity along the x direction, it is suitable to 
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express x and x as 

x = nC, + s, x = n £ + s , (13) 

with n,n £ Z and s, s £ [0, £). Since a(x) is peri- 
odic with wavelength £ it follows that a(x) = a(s) and 
x — x = (n — n )C + s — s depends on n and n only 
via the difference n — n. This property and translational 
invariance along the y direction imply the Fourier decom- 
position 

OO p p 

M(p x ,p y ;s,s ) = ^2 d V d y 

n,n — — oo 

x M(n -n\y-y';s, s ) e - i P^ n - n ')« e - ip ^ y - y ^ (14) 

where M(p x ,p y ; s, s ) is diagonal. Furthermore with s, 
s G [0, C) one can form 

C k i {PxiPy) = [ ds f ds 
Jo Jo 

x e"^M( Pl ,p 9 i S ,;)e- Ms ' /( . (15) 



As expected, the different representations of M in 
terms of the matrices C or B (C is the Fourier transform 

of M and B is the Fourier transform of e" ,!,lS Me ,PlS ) do 
not change the final result for the force [Eq. IjlUI) ] which 
is given by 0] 




xtr(B- 1 {p x ,p y )dB{ Px ,p y )). (16) 



Here £r denotes the partial trace with respect to the in- 
dices k, I of the infinite-dimensional matrix B (or C) and 
we have taken into account the contribution of both fluc- 
tuating components of the director field. 

In the following we continue with the block-diagonal 
form of the matrix kernel M since in this representation 
the patterning functions and bk are somehow simpler 
than their counterparts in the diagonal form of M. Ac- 
cordingly, the matrices N m [Eq. l(T^)l ] are given by 



N m = 



/C(o) K b M 

\<t%W $*(o) 

A& — A a 
2A Q 



Sm,0 



2p 



KP (l-At,p) J 
2 2p 



(1-At,p) 2 -pd zg Agp 
2p C 2p 2 



a (f -A 6 )+6 (f -A 6 ) 



«o(f -A b )+6o(f -A b ) 



- pd 



for m even, and 

\ — A fl 



2A a 



-A b ) e -f rf + a m (^ 



A h e"P" 



- \ b )e~P d + b m (^~ A fc )e-P- 



(17) 



(18) 



for m odd, with z a(6) = 1 + A " Aa A " ao(Mi ^(fc) = 

1 + ^IT^MM, Pm = ^(p x + 2Trm/0 2 +pl, P = 
p m =o, a = b = 1/2, a m=odd = (-l)(l^l-l)/2/( 7r | m |) j 
b m = odd = e- 2 ™" lA /<a m , and 

«'w=(^) 2 f:v™ do) 



where the prime at the summation sign indicates that in 
the sum the terms with even k are excluded. 



III. LATERAL FORCES 



A. Fluctuation-induced shear forces 



In the limit d 3> £ we find that the contributions from 
the elements Bki [Eqs. (|llf> and (|12[) ] to the force de- 
crease rapidly with increasing absolute values of fc, Z, so 
that the expression for the force [Eq. I|l(j|) ] converges al- 
ready at small orders of M with k, I — —M, • • • , M. 
Taking into account only the elements Bki [Eq. I|ll|l ] 
for k,l = —1,0,1, the asymptotic behavior of the 
fluctuation-induced lateral force J^ut = —dgF in the limit 
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| At - A Q |/A a < 1 is given by 

•fiat(rf>C) _ 8(A fc - A g ) 2 27rd/c /2tt5n 



x/(d/C,A b /C) + 



with the dimcnsionless function 
f(u,v) = 



v C / 



A/, — A n 



(20) 



(1 + 2™) 2 Jo 



i l+2 t ,(l-t»^)^+(^) 2 (l+47rt.) ^ 2ti3 . 

1 (T^P e 



1 



/ 1+vx \ ' 
\ 1—vx J 



(21) 



The lateral force oscillates as function of 8 reflecting 
the underlying lateral periodic pattern and its magni- 
tude decays exponentially as function of d/C, because 
/(oo, v) < CO. 

For arbitrary values of d, we evaluate the force in 
Eq. I|16fl numerically. Although the matrix B [Eq. (|ll(l ] 
is infinite-dimensional, the value of the force saturates 
at some finite values for k,l. Our numerical results for 
the fluctuation-induced lateral force as function of the 
shift 8 for arbitrary strength of the contrast A;, — A Q 
are shown in Fig. [21 The lateral force acts against 
the increase of the lateral displacement 8 in the inter- 
val [0, C/2) by being a restoring force and acts favorably 
with the increase of 8 in the interval (C/2, C] by being 
a pulling force. Therefore the force is antisymmetric 
with respect to §/(, — 0.5. Upon approaching the max- 
imum misalignment, i.e., 8/C, = 0.5, the restoring force 
vanishes. This implies that the interaction free energy 
Vi&t(8) — ~ J 3~\ a t(S )dd has its maximum at S/C = 0.5 
where the opposing parts of the substrates face each other 
and attains its minimum at 8 = 0. The force is maximal 
at S/C = 0.25 and 0.75 where Mat exhibits its strongest 
dependence on 8 (Fig. |3J. In Fig. 0] we show the decay 
of as function of d. Asymptotically, T\ at /(^sTS'/C 3 ) 
vanishes as -0.003 e - 27r!i /<7(d/C) for A Q /C = 4, A fc /C = 8, 
and S/C = 0.25. 

We note that due to the assumption £ a = £fc inter- 
changing A a and Ah leaves the system unchanged [Fig. 2]. 
Thus the force must be identical for A a «-» Ah. While in 
Eqs. I|20|) and l|21l) this symmetry is explicitly valid up 
to the second order in (Ah — A a ), the numerical results 
respect this symmetry fully. This provides a very useful 
check of the numerical calculations. 




d/C = 0.01 -m- 

d/C=0.02 
d/C = 0.05 -•- 
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0.75 



FIG. 2: The fluctuation-induced lateral force ^i at in units 
of ksTS/^ 3 between two periodically patterned substrates at 
distance d as function of the shift S in units of the periodicity 
C for d/C = 0.01, 0.02, and 0.05 (see Fig. [TJ. The anchoring 
on the stripes in terms of the extrapolation lengths is taken 
to be X a /C = 4 and Aj/C = 8. J-ut is antisymmetric around 
5/C, = 0.5. There is no restoring force if the misalignment is 
maximal, i.e., at S/( = 0.5. J~ut < means that the plates 
are pulled back towards the preferred alignment at 8 — 0; 
for T > the plates are pulled forward towards preferred 
alignment at 8 = £. 
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FIG. 3: The effective lateral potential VJat(<5) = 
— Jq Jiat(<5 )d8 in units of kBTS/C, 2 between two periodically 
patterned substrates at distance d as function of the shift 8 in 
units of the periodicity £ for d/C, = 0.02, 0.03, and 0.05 and 
X a /C — 4, Xt/( = 8. For all values of d/C, the inflection points 
of the potential are at 8/C, = 0.25 and 0.75. 



B. Lateral van der Waals force 

Endowing the substrates with the envisaged stripe pat- 
terns requires corresponding chemical patterns which in 



turn involve at least two different species providing the 
chemical contrast. These species do not only interact 
(differently) with the nematic liquid crystal in the vicin- 
ity of the substrate, giving rise to two different extrapola- 
tion lengths A a and Ah, but also interact across the liquid 
crystal with each other via dispersion forces. The later 
interaction provides a lateral force as well which adds to 
the fluctuation-induced lateral force. Note that such a 
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ters. This leads to the following expression for the van 
der Waals potential energy between the two monolayers: 



vdW 



Iat 



with 



and 



d 2 x, 



f ,2 E(x 1 ,x 2 ;S) 
Is X2 [d 2 + ( Xl -x 2 ) 2 ]3 



£7(2:1,0:2; S) = E AA a{xi)b{x 2 \8) + 
E BA [a(xij - 2a(x 1 )b(x 2 ;S) + b(x 2 ; 6)] 
+E BB [l-a{x 1 )][l-b{x 2 -8)] 



(23) 



(24) 



FIG. 4: The fluctuation-induced lateral force .Fiat in units 
of UbTS/C, 3 between two periodically patterned substrates 
as function of the film thickness d in units of the period- 
icity (. The anchoring on the stripes in terms of the ex- 
trapolation lengths is taken to be \ a /C = 4, Xb/C — 8, 
and the lateral shift between the patterns is d/C, = 0.25. 
Asymptotically, Fi at /(fcsT5'/^ 3 ) for d/C, = 0.25 vanishes as 
-0.003 e" 27rd/c /(d/C) (dashed line). The variation of F Jat 
for small values of d/C, is shown in the inset for A a /C — 4 
(circles) and X a /C — 2 (triangles). The other system parame- 
ters remain the same. As expected the absolute value of the 
amplitude of the lateral force increases upon increasing the 
contrast I At. — AJ. 



lateral force due to direct interactions is the same if the 
two patterned substrates are separated by vacuum or by 
a nematic liquid crystal, as long as the mean nematic 
order is not affected by the stripe patterns. 

In order to estimate these direct interactions between 
the patterned substrates, we consider each substrate to 
be covered by a monolayer whose chemical composi- 
tion varies periodically, alternating between A- and B- 
particles. We neglect non-additivity aspects of the dis- 
persion forces and consider pairwise interactions between 
the patterned monolayers at z = and z — d. Since we 
consider d to be large compared with the diameters of 
the A and B particles, we can disregard that particles 
forming the monolayers occupy discrete lattice sites. As 
pair potentials between the two species we take Lennard- 
Jones potentials 



Uij (r)=4e ij [(^) -(2L)],i,j=A,B (22) 

where A(B)-particles give rise to the extrapolation length 
A a (6) . Since in the present context d ^> ov, , for the lateral 
force only the attractive part of the pair potentials mat- 



Eij - Ae i ja ij 'E i 'E j 



(25) 



where Tia(b) is the areal number density of A(B)- 
particles in the monolayer forming the stripe A a (;,). Car- 
rying out the integration over ?/i and y 2 in Eq. (|23l) . one 
obtains 



V! 



vdW 
lat 



(5) = / dxi [ dx 2 

° J-L/2 J-L/2 

E(xi,x 2 ;S) 

x [d^ + { Xl ~x 2 fY>/^ 



(26) 



where L is the lateral extension of the system both in the 
x and y directions. From this the lateral van der Waals 
= -~ d sVilt W ( S ) can be calculated: 



force 



rvdW 
•Mat 

S 



E 



r w (5/c,d/o 



(27) 



with the scaling function 



_^ r n+l/2+S/C 

s/C 



rl/2 

< / dX 2 
10 



X 2 — X\ 



[(d/C) 2 + (X 2 - Xi)2]V 



dXi 

2 (M) 



where E = Eaa — 2Eab + Ebb- The force and the 
potential as function of 8/C, are shown in Fig. |SJ The 
comparison between Figs. |21 and |S] (a) reveals that the 
lateral van der Waals force is practically constant over 
a wide range of shift values and varies steeply around 
the positions of maximum and minimum misalignment 
while the fluctuation-induced force varies more smoothly 
across all shift values. Apart from that, the qualitative 
features are the same for both forces. 



In order to estimate the lateral van der Waals force 



we assume that the particles are closely packed within 
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0.25 0.5 0.75 1 0.25 0.5 0.75 



s/C s/c 

FIG. 5: (a) The dimensionless scaling function f vdw [Eqs. 1271 and 1281 1 of the lateral force induced by direct van der 
Waals interactions between the chemically patterned monolayers covering the two substrates at distance d as function of the 
shift S in units of the periodicity C of the pattern for d/C = 0.05 and 0.06. (b) Corresponding lateral potential V{^ W (S) = 



^latjj ~ 3 P N S/(Mm) 2 for the contribution cx E^ 
actual lateral van der Waals force JYat 



the monolayer |l7j| . so that the areal number density 
is E A(B) = %/3/(6i? 2 ) w 2[<t aa(bb) ]- 2 /V3 where R = 
gaaI^ {obbI'2) is the radius of the A [B) particles form- 
ing a triangular lattice. According to Table I in Ref. |is[ 
typical values are e^/fes ~ 300 K and cr^ « 0.5 nm. For 
C = 200 nm one has Eij/C 5 ~ L5 x 10 ~ 5 pN/(^m) 2 . 
According to Fig. (a) for d/C = 0.05 this implies 

to the 
proportional to 
E = Eaa — 2Eab + Ebb- Thus for a suitably chosen 
constant E, without compromising the goal of achieving 
A a = Aft/2 (as used in our calculations), the lateral van 
der Waals force can be quite smaller than 3 pN per area 

s/(H 2 

From Fig. [21 one finds, for the same system param- 
eters C an d d considered above and at room temper- 
ature, for the fluctuation-induced lateral force J-"i a t ~ 
0.02 pNS/(/xm) 2 . Thus the background lateral van der 
Waals force tends to be stronger than the nematic fluc- 
tuation induced force. However, for a suitably chosen 
chemical contrast of the particles forming the chemi- 
cal stripes, it appears to be possible to determine the 
fluctuation-induced lateral force by measuring the shear 
force once with and once without the nematic liquid 
between the patterned substrates. The ratio of the 
two forces for d/C = 0.05, A Q /C = 4, \ b /( = 8, 
and E/C b w 1.5 x 10~ 5 pN/(/im) 2 as function of S 
is shown in Fig. It appears that the fluctuation- 
induced lateral force becomes more prominent around 
S/C = 0.25,0.75. From Figs. |3 and (a) one no- 
tices that both J-j a t/S and J 7 ^™ / S vanish linearly at 
S/C = 0,0.5,1.0 [as (-0.084 + 0.168 S/C) pN/(//m) 2 and 
(-33.5 + 66 S/C) pN/(^m) 2 at S/C = 0.5, respectively, for 
d/C = 0.05 , E/C 5 = 1.5 x 10~ 5 pN/(^m) 2 , £ = 200 nm, 
and k B T = 4 x 10~ 21 J] and since the slope of J^ dw 



lat 



is much larger than the corresponding slope of J-"i a t, the 
ratio ^i at /^ w at S/C = 0, 0.5, 1.0 is small. From Fig.© 
one should not draw the conclusion that the lateral fluc- 
tuation induced force is at most half a percent of the cor- 



0.005 
0.004 

0.003 

(j^ 0.002 
0.001 





• 1 • 







25 



0.5 

s/C 



0.75 



FIG. 6: The ratio of the fluctuation-induced lateral force _7-i at 
and the lateral van der Waals force J-JJIt as function of the 
shift S in units of the periodicity ( for d/C = 0.05, A /C = 4, 
\ b /C = 8, E/C 5 « 1.5 x 10~ 5 pN/(/im) 2 , and at T = 290 K. 
T\ax is more prominent around S/C = 0.25 and 0.75. 



responding lateral van der Waals background force. This 
ratio is inversely proportional to E = Eaa—2Eab+Ebb- 
Figure El corresponds to a parameter choice for which E 
is estimated by an individual Eij and not by the actual 
contrast expressed by E, which vanishes for A = B. Ac- 
cordingly, for suitable choices of A and B, E can be sig- 
nificantly smaller than the Eij used in Fig. which then 
leads to a significantly larger ratio. 



IV. FLUCTUATION-INDUCED NORMAL 
FORCE 

The effect of a periodic anchoring at one substrate, 
with the second substrate being homogeneous, on the 
fluctuation-induced normal force was studied in Ref. |ll| . 
It turned out that for the description of the normal force 
the single patterned substrate can be replaced by a uni- 
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d/c s/( 



FIG. 7: The fluctuation-induced normal force f mrln in units 
of knTS/C, 3 between two periodically patterned substrates as 
function of the film thickness d in units of the periodicity C- 
The anchoring on the stripes in terms of the extrapolation 
lengths is taken to be Xb/C — 8 and X a /C = 4 (full line), 2 
(dashed line). For d < X a ,Xt as shown here the anchoring 
is weak but finite at both substrates and the force is at- 
tractive. As expected the absolute value of the amplitude of 
the force decreases upon decreasing the extrapolation length. 
The force depends very weakly on the shift 8. Here 8/C is set 
to 0.3. 

form substrate with an effective anchoring strength, i.e., 
the force is given by the force found between two uniform 
substrates characterized by their effective anchoring. De- 
pending on the model parameters, the normal force is ei- 
ther repulsive or attractive - corresponding to an effective 
similar-dissimilar or an effective similar-similar boundary 
condition, respectively. In the present case of two pat- 
terned substrates, we have calculated the normal force 
numerically [Eq. i|16|) ]. Figure [7| shows the dependence 
of the normal force f nonn on d. For d < X a , Ab as shown 
here, the force is attractive and decays monotonically as 
function of d. In this regime, the anchoring is weak at 
both substrates so that the effect of the periodicity is 
not visible. We note that in this case the fluctuation- 
induced normal force [Fig. [7] is about 1000 (100) times 
larger than the fluctuation-induced lateral force [Fig. 0] 
for X a /C — 4 (2). Figure however, shows the behavior 
of ^-"norm as function of 6 in the regime of strong but finite 
anchoring d > A a , A/, . In this case the normal force is os- 
cillatory and attractive, and its magnitude is comparable 
with the fluctuation-induced lateral force [Fig.|5]. 



FIG. 8: The fluctuation-induced normal force f mrm in units 
of feflTS/^ 3 between two periodically patterned substrates at 
distance d/C = 0.45 as function of the shift 8 in units of the 
periodicity f . The anchoring on the stripes in terms of the ex- 
trapolation lengths is taken to be X a /C = 0.1 and Xb/C = 0.3 
(triangles), 0.4 (circles). For d > X a ,Xb as shown here the 
anchoring is strong but finite at both substrates and the force 
is attractive. As expected the absolute value of the ampli- 
tude of the force decreases upon increasing the extrapolation 
length The force oscillates as function of the shift 8 and 
for complete misalignment, i.e., 8/C = 0.5 the attraction is 
weakest . 




0.1 0.15 0.2 0.25 0.3 0.35 0.4 



(h - K)/C 

FIG. 9: Fluctuation-induced normal force jT norm divided by 
fluctuation-induced lateral force J-i^t as function of the con- 
trast At — A a in units of the periodicity C f° r X a /C = 0.1, 
d/C = 0.45, corresponding to strong but finite anchoring, and 
8/C = 0.25. J--norm and ^i at are of comparable size. 



V. SUMMARY AND CONCLUSION 

We have calculated the fluctuation-induced forces act- 
ing on two substrates chemically modulated with period 



£ and confining a ncmatic film of thickness d (Fig. 
The substrates are characterized by homeotropic anchor- 
ing with alternating extrapolation lengths A a and Afc. We 
have studied the shear force as function of the lateral 
shift S between the patterns on the substrates and of 
their separation d. For |A& — A a | <C A a and £ <C d, 
the lateral force sinusoidally oscillates as a function of 
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5/C and decays exponentially with d/Q [Eq. I|20|l and 
Fig. 0]. For stronger contrasts, the lateral force and its 
corresponding potential have been evaluated numerically 
[Figs. [21 and [3]. It turns out that for a suitably chosen 
chemical contrast the fluctuation-induced lateral force is 
comparable [Fig. E] with the background lateral van der 
Waals force between the corresponding monolayers on 
the top and bottom substrates forming the chemical het- 
erogeneity [Fig. [5] . In order to complete the picture of 
the forces in the presence of the patterned substrates we 
have also calculated numerically the fluctuation-induced 
normal forces [Figs. [7| and |H] and found them to be com- 
parable with the fluctuation-induced lateral force in the 
case of strong anchoring [Fig. E] . 

Patterning at small length scales gives rise to rich in- 
terfacial phenomena. Surface modulations open the pos- 
sibility of controlling the morphology of wetting films and 
generate structural phase transitions which are central to 
the behavior of structural forces induced by distortions of 
the liquid crystal order parameter. In such cases, in addi- 
tion to the fluctuation-induced forces, the substrates are 
subject to liquid-crystalline elastic forces Q. Such elastic 
forces scale with K [Eq. @] and are for large d larger than 
the fluctuation-induced forces which scale with fceT. 
Since for small d the elastic forces scale as d~ 5 |2| but 



the fluctuation-induced force as d~ 3 p], the latter can, 
however, even dominate for small d. The mean field di- 
rector contributions are completely eliminated for such 
model parameters and boundary conditions for which the 
director structure is uniform. In Ref. [19j it has been 
demonstrated that this uniformity can indeed occur for 
suitable combinations of model parameters, even in cases 
of competing planar and homeotropic anchoring condi- 
tions, which have not been considered here. Capillary 
forces due to capillary condensation and the formation 
of bridge phases [2(], 0, are other sources for the 
structural forces in the vicinity of the nematic-isotropic 
phase transition. If the fluid is confined to very narrow 
slits, patterning may give rise to capillary bridges of dif- 
ferent liquid crystalline order. Under such conditions, it 
would be interesting to study the stress under the shear 
strains by shifting the lateral substrate structures out of 
phase 23] which might give rise to rather strong lateral 
forces. 

Lateral forces can technologically be used to align the 
parallel substrate structures. For instance, for those 
ranges of the model parameters for which the liquid crys- 
talline lateral forces are significant, the liquid crystal can 
be filled into the slit pore to align the substrate structure 
and then be removed. 
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